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TRANSCUTANEOUS DETECTION OF ARTERIAL BLEEDING BY TWO- 
DIMENSIONAL ULTRASONIC IMAGING OF TISSUE VIBRATIONS 

Government Rights 

This invention was funded at least in part with grants from the U.S. Office of 
5 Naval Research (ONR) (Award Nos.N00014-96-l-0630, N00014-01-G-0460, and 
N00014-99-1-0982), and grants from the U.S. Department of Defense (DoD) (ARMY 
MRMC Award Nos. DAMD 17-02-2-00 14, and DAMD 17-00-2-0063), and the U.S. 
government may have certain rights in this invention. 

Field of the Invention 

10 This invention generally relates to a method for identifying tissue vibrations 

using an ensemble of two-dimensional (2D) ultrasound data acquired with an ultrasound 
transducer, and more specifically, to a method and system for identifying bleeding sites 
within a body of a patient in real time, based upon the tissue vibration that occurs at the 
sites, so that arterial bleeding and bleeding from organs can readily be localized in an 

15 image. 

Background of the Invention 

Internal bleeding is a significant cause of death in cases of trauma, and rapid and 
effective diagnosis of patients with uncontrolled bleeding has long been recognized as an 
important goal to lower mortality and morbidity. Currently, internal bleeding is 

20 diagnosed using angiography, Computed Tomography (CT), diagnostic peritoneal 
lavage, B-mode ultrasound, and exploratory laparotomy, while surgical intervention is 
the most common treatment option. The detection of a bleeding site (i.e., diagnostic 
imaging) and the closure of a bleeding wound should preferably be accomplished as 
quickly as possible to minimize blood loss in a patient and reduce mortality and 

25 morbidity associated with such blood loss. 

Advances in duplex and color-flow ultrasound in the last two decades have had a 
significant clinical impact on vascular diagnosis. For example, the use of Doppler 
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ultrasound has been shown to be effective for targeting a bleeding site, as disclosed in a 
paper by R.W. Martin, S. Vaezy, P. Kaczkowski, G. Keilman, S. Carter, M. Caps, K.W. 
Beach, MI. Plett, and LA. Crum, entitled "Hemostasis of punctured vessels using 
Doppler-guided high-intensity ultrasound," Ultrasound Med Biol, vol. 25, pp. 985-990, 

5 1999. However, this technique suffers from the disadvantage of imaging a limited region 
of interest. Although color-flow ultrasound can image a large region of interest, currently 
it lacks sufficient sensitivity for diagnosing internal bleeding, due to weak scattering from 
blood and the slow flow velocity of blood bleeding from a wound, especially in the case 
of deep bleeds and organ bleeds. In a paper by X. Shi, R.W. Martin, S. Vaezy, and L.A. 

10 Crum, entitled "Quantitative investigation of acoustic streaming in blood," X Acoust 
Soc. Am., vol. Ill, pp. 1110-1121, 2002, the use of acoustic streaming is proposed for 
distinguishing between stagnant blood and tissue using color-flow images. However this 
paper does not suggest how to detect a bleeding site with acoustic streaming. The use of 
contrast agents has also shown to be promising for localizing active bleeding sites, as 

15 discussed in a paper by J-B Liu, D.A. Merton, B.B. Goldberg, N.M. Rawool, W.T. Shi, 
and R Forsberg, entitled "Contrast-enhanced two- and three-dimensional sonography for 
revaluation of intra-abdominal hemorrhage," J. Ultrasound Med, vol. 21, 
pp. 161-169, 2002. Yet, the use of contrast agents is time-consuming and can sometimes 
be dangerous. 

20 Each of these prior art techniques for detecting bleeding is unable to provide 

efficient real-time images in which the location of a bleeding site can be rapidly 
identified. Also, as discussed in greater detail below, simply using conventional 
color-flow data for imaging a site cannot readily distinguish between pooled blood and 
bleeding at the site. Accordingly, it is necessary to develop a new method for rapid 

25 diagnosis of internal bleeding. 

Physical examination is an important element of the initial assessment of a 
trauma patient with suspected internal bleeding. If an audible "bruit" or a palpable 
"thrill" is found upon physical examination, further diagnostic tests for internal bleeding 
or surgical intervention are often recommended. It has now been established that bruits 

30 and thrills are produced by the forces exerted on vessel walls by eddies created as blood 
flows from a high-pressure region to a low-pressure region through a narrow orifice. The 
pressure fluctuations in eddies cause local vibrations in the vessel wall and surrounding 
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tissue and manifest either as bruits or thrills at the skin surface. The power spectrum of 
the vibrations exhibits a frequency peak called the "break frequency," which is directly 
related to the diameter of the orifice and the local flow velocity through the Strouhal 
number. In conventional color-flow ultrasound images, tissue vibrations from abnormal 

5 blood flow produce characteristic speckled artifacts in the surrounding tissue. But, these 
artifacts are difficult to interpret and are not quantitative. Tissue vibrations have been 
previously studied using one-dimensional (ID) pulsed Doppler ultrasound, and the prior 
art includes disclosure of a wavelet-based method for detecting and characterizing arterial 
vibrations in internal bleeding (M.I. Plett, "Ultrasonic arterial vibrometry with 

10 wavelet-based detection and estimation," PhD. dissertation, Univ. of Washington, 2000). 
However, this pulsed Doppler-based technique also has a limited field of view, and is 
along a single scan line. Furthermore, as disclosed in this paper, the processing was done 
off-line, so it was not possible to create images interactively in real time. 

Accordingly, it would be desirable to develop a new tissue vibration detection 

15 and imaging mode for ultrasound instruments in which vibrations produced by blood 
flow eddies can be detected and color-coded according to their amplitude and frequency 
and overlaid on a B-mode and/or a color-flow image in real time. The tissue vibration 
imaging mode might then be used for locating the origin of the vibration more precisely, 
relative to the patient's anatomy and/or for obtaining simultaneous information about 

20 vibrations and the underlying blood flow. Acoustic hemostasis using High Intensity 
Focused Ultrasound (HIFU) is a promising new technique for stopping internal bleeding 
without invasive surgical intervention. For effective targeting and monitoring of 
hemostasis, non-invasive real-time localization of a bleeding site in real time is essential. 
This new technique might thus be used both for diagnostic determination of a bleeding 

25 site, and optionally, in conjunction with HEFU or other desired therapy, for localizing the 
bleeding site in real time, so that therapy might be effectively applied to stop the bleeding 
as indicated in the U.S. Patent No. 6,083,159, "Methods and devices for providing 
acoustic hemostasis", U.S. Patent No. 5,993,389, 'Devices for providing acoustic 
hemostasis", and US patent 5,882,302, "Methods and devices for providing acoustic 

30 hemostasis". 
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Summary of the Invention 

This invention was developed to localize a bleeding site by imaging tissue 
vibrations caused by arterial bleeding into a free space. Sounds caused by these 
vibrations (bruits) are sometimes audible using a stethoscope or palpable at the skin 

5 surface and are indicative of internal bleeding in a trauma patient. Previously, a 
Doppler-based algorithm might have been used for offline analysis of vibrations along a 
single ultrasound scan line. . In contrast, this invention employs an algorithm that 
processes an ensemble of 2D ultrasound data for detecting and imaging tissue vibrations 
in a relatively large region of interest. The algorithm has been implemented in a 

10 programmable ultrasound system to study the usefulness of tissue vibrations in real-time 
localization of bleeding sites in animals and has been shown to be very effective for this 
purpose. 

The vibration imaging algorithm uses an ensemble of 2D ultrasound data 
acquired during conventional ultrasonic imaging and utilizes the clutter signal, which is 

1 5 normally suppressed in conventional color-flow imaging, to detect and characterize tissue 
vibrations. Three signal processing algorithms have been developed that are suitable for 
this purpose, including one based on phase decomposition, a second based on spectral 
estimation using eigendecomposition, and a third based on spectral estimation using 
autoregressive modeling for isolating vibrations from clutter, blood flow, and noise in an 

20 ensemble of 2D ultrasound data. 

Real-time tissue vibration imaging has been implemented at frame rates, for 
example, of 10 frames/second, on an ultrasound system with a software-programmable 
signal and image processing back-end. The preliminary results confirm that vibrations 
produced as a result of arterial bleeding can be detected and imaged with this invention. 

25 The vibration amplitude is expected to be the largest near the site of the bleeding, and this 
fact can be used to localize a bleeding site quickly and non-invasively. The strong 
backscattered ultrasonic echoes from tissue vibrations can improve visualization of 
internal bleeding sites that are otherwise hard to image due to weak scattering from 
blood. 

30 Potentially, this new tissue vibration imaging technology could be useful in a 

variety of devices and clinical settings. For example, a low-cost portable screening 
device with tissue vibration detection functionality could be beneficially employed by 
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paramedics and trauma centers to localize bleeding in patients. In addition, a tissue 
vibration imaging mode on high-end ultrasound systems can augment duplex ultrasound 
for enhanced diagnostic capability. The detected tissue vibrations indicative of internal 
bleeding could be presented as an audible signal in a manner recognizable to a person 
5 trained to listen to bruits using a stethoscope or as a palpable signal recognizable to a 
person trained to detect palpable thrills. 

One aspect of the invention is directed to a method for detecting and localizing 
internal bleeding using an ensemble of 2D ultrasound data by detecting and 
characterizing tissue vibrations caused by blood flow eddies at the bleeding site. The 
10 method includes the step of processing an ensemble of 2D ultrasound data, producing a 
tissue motion spectrum signal at a site. The tissue motion spectrum signal is then filtered 
to produce a filtered signal from which any contribution to the tissue motion from a 
source other than bleeding at the site has been substantially minimized, producing a 
tissue vibration signal. A vibration image is displayed using the filtered signal and 
1 5 indicates a location of bleeding at the site. 

One approach for processing the ensemble of 2D ultrasound data comprises the 
steps of estimating a correlation matrix from the ultrasound data, and carrying out an 
eigendecomposition of the correlation matrix to identify a signal subspace and a noise 
subspace. A frequency of the dominant vibration components in the signal subspace and 
20 the noise subspace are then estimated, and based upon that estimate, a vibration 
amplitude estimate and a vibration frequency , estimate are determined. At least one of 
the vibration amplitude estimate and the vibration frequency estimate comprises the 
tissue vibration signal. 

A second approach for processing the ensemble of 2D ultrasound data includes 
25 the step of computing reflection coefficients from the ultrasound data. Linear prediction 
filter coefficients are computed from the reflection coefficients. A power spectrum is 
estimated, and the peaks in the power spectrum are detected. Based upon the estimate of 
the power spectrum and the peaks, a vibration amplitude estimate and a vibration 
frequency estimate are determined. At least one of the vibration amplitude estimate and 
30 the vibration frequency estimate again comprises the tissue vibration signal. 

In yet another approach for processing the ensemble of 2D ultrasound data, a 
mean clutter velocity is estimated from the ultrasound data using autocorrelation. The 
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ensemble of 2D ultrasound data are down-mixed with the mean clutter velocity, 
producing a down-mixed signal. A phase of the down-mixed signal and a mean phase of 
the down-mixed signal are determined, and the mean phase is subtracted from the phase 
of the down-mixed signal, producing a residual phase. Hie residual phase is then 
5 decomposed into its dominant components. By applying energy and frequency 
thresholds, any contribution to the tissue vibration due to noise and blood flow are 
substantially suppressed, yielding an estimate of vibration amplitude and vibration 
frequency of tissue due to bleeding. 

The step of decomposing the residual phase preferably comprises the steps of 
10 estimating a correlation matrix from the residual phase, and performing an 
eigendecomposition of the correlation matrix to determine the dominant components. 

The step of filtering preferably comprises the step of filtering out clutter and 
noise at frequencies that are substantially lower than an expected frequency range of 
tissue vibrations corresponding to bleeding at the site, and also preferably includes the 
15 step of filtering out noise that is at frequencies, which are substantially higher than an 
expected frequency range of tissue vibrations corresponding to bleeding at the site. 

In addition, the method can include the step of confirming that vibrations 
displayed in the vibration image correspond to bleeding at the site, by placing a Doppler 
sample volume at a location of the tissue vibration determined from the tissue vibration 
20 image. In this step, the tissue vibration spectrum determined from a relatively larger 
number of samples (typically, 64-512) available using Doppler data can be employed to 
confirm the tissue vibration detected from only a few ensemble (typically, 6-16) of 2D 
ultrasound data. 

The step of displaying the vibration image preferably comprises the step of 
25 displaying at least one of a vibration amplitude image and a vibration frequency image of 
the site. Because of its efficiency, the method can include the step of displaying the 
vibration image in connection with an underlying anatomy of the site (Le., the B-mode 
grayscale image), substantially in real time. 

Another aspect of the present invention is directed to apparatus for detecting and 
30 localizing bleeding at an internal site using an ensemble of 2D ultrasound data. The 
apparatus includes an ultrasound transducer for transmitting ultrasound pulses toward the 
internal site and receiving ultrasound data from scatterers at the internal site, including 
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tissue that is vibrating due to bleeding at the site. In one embodiment, a front-end system 
controls the ultrasound pulses produced by the ultrasound transducer and demodulates 
the echoes received by the ultrasound transducer, producing a color-flow signal having 
both in-phase and quadrature components. The apparatus also includes a back-end 

5 system to receive the color-flow signal from the front-end system and which includes a 
tissue vibration processor. The front and back-end systems can optionally be combined 
into a single unit, or one or more parts of these systems can be operating remotely from 
other parts of these systems. The tissue vibration processor processes the ultrasound 
signal to estimate tissue vibrations caused by bleeding, producing a tissue vibration 

10 signal. In one embodiment, the tissue vibration signal is converted to an image signal by 
the back-end system. A display is coupled to the back-end system to receive the image 
signal, to display a tissue vibration image in which bleeding at the internal site is 
indicated Optionally, the display could be decoupled from the tissue vibration detection 
and identification apparatus. For example, the display can be located physically in a 

15 hospital, while the tissue vibration detection and identification apparatus is located 
physically at another location, such as in an ambulance. Also, instead of a visible 
display, the result of tissue vibration detection and identification can be presented as an 
audible or a palpable output indicating tissue vibrations. The tissue vibration signal can 
also be interpreted by an automated algorithm to indicate internal bleeding, and the result 

20 of the automated interpretation presented as an electronic readout Generally, the 
functions performed by the apparatus are consistent with the steps of the method 
described above. 

Brief Description of the Drawing Figures 

The foregoing aspects and many of the attendant advantages of this invention will 
25 become more readily appreciated as the same becomes better understood by reference to 
the following detailed description, when taken in conjunction with the accompanying 
drawings, wherein: 

FIGURE 1 is a functional block diagram of an exemplary ultrasound system that 
is suitable for carrying out tissue vibration imaging in accord with the present invention; 
30 FIGURE 2 is a graph of an expected Doppler spectrum (i.e., frequency vs. 

power) from a vibrating sample volume; 
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FIGURE 3A is a flow chart showing the logical steps of a subspace-based 
algorithm for creating a vibrating tissue image in which bleeding is evident; 

FIGURE 3B is a flow chart showing the logical steps of an alternative algorithm 
that uses autoregression for creating a vibrating tissue image in which bleeding is 
5 evident; 

FIGURE 3C is a flow chart showing the logical steps of yet another alternative 
algorithm that uses phase decomposition for creating a vibrating tissue image in which 
bleeding is evident; 

FIGURE 4 is a schematic diagram showing how a gel standoff is used for 
10 acoustically coupling an ultrasound transducer to an organ that is exposed and incised to 
create a bleeding site; 

FIGURES 5A, 5B, and 5C are respectively a color-flow image of a punctured 
femoral artery, a vibration amplitude image of the artery, and a vibration frequency 
image of the artery; 

1 5 FIGURES 6A, 6B, and 6C respectively illustrate a color-flow image of an incised 

spleen, a vibration amplitude image of the spleen, and a vibration frequency image of the 
spleen; and 

FIGURES 7 A, 7B, and 7C respectively illustrate a color-flow image of an incised 
liver, a vibration amplitude image of the liver, and a vibration frequency image of the 
20 organ. 

Description of the Preferred Embodiment 

Tissue Vibration Imaging System 

FIGURE 1 is a block diagram illustrating an ultrasound system 10, which is 
generally similar to a conventional ultrasound system, but which has been modified to 

25 include tissue vibration imaging and is thus usable in practicing the present invention. 
Ultrasound system 10 includes an ultrasound transducer 12 that transmits a signal, which 
is modulated with a carrier frequency, typically 1 MHz - 15 MHz, using multiple cycles 
(i.e., 2-20 cycles). The transmitted signal is reflected by scatterers (not shown) along the 
beam path and is received after a time delay, which depends upon the distance of the 

3 0 scatter from the transducer. In an acquisition stage, the acoustic echoes received from the 
tissue are converted to electrical signals by the transducer, and these signals are digitized 
by analog-to-digital converters (not separately shown). A front-end subsystem 14 
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includes a beam former 16 that performs dynamic focusing, apodization, and steering of 
both transmitted and received ultrasonic waveforms. Also included in front-end 
system 14 are a time-gain-compensation (TGC) circuit 18 that amplifies signals with a 
variable gain that is proportional to the depth within tissue, and a radio frequency (RF) 
5 demodulator and decimator20 that digitally removes the higk frequency carrier by 
quadrature demodulation and decimation, providing both in-phase (I) and quadrature (Q) 
samples, which may be represented as a complex quantity, I(t) + jQ(t). The acquired 
quadrature ensemble (or color-flow) data are then processed in a back-end subsystem 22, 
depending on the one (or more) ultrasound mode(s) that is/are selected, e.g., B-mode, 
1 0 color-flow mode, tissue vibration mode, and Doppler mode. 

For producing anatomic images of tissue, the signal of interest is the envelope of 
I(t) + JQ(t). A B-mode processor 24 computes the magnitude of the echo, 
B a (t) = ^I 2 (t)+Q 2 (t) and compresses the dynamic range to make it suitable for display as 
a grayscale image on a monitor. The time delay introduced by the scatterers is reflected 
15 in the phase of the complex quantity I(t) + jQ(t). Thus, the phase of the complex 
received signal provides an estimate of the instantaneous position of the scatterer. By 
monitoring the change of phase over time, the displacement and velocity of the scatterer 
can be estimated. In color-flow imaging, multiple pulses (commonly from 6 to 
16 pulses) are transmitted and received along each scan line at a rate known as the pulse 
20 repetition frequency (PRF). A collection of received temporal samples from each spatial 
location is thus called an "ensemble." A color-flow processor 26 estimates the blood 
flow velocity from the ensemble of data by estimating the phase difference between the 
adjacent temporal samples, typically using an autocorrelation algorithm. A 2D image is 
created by acquiring multiple samples from different spatial locations. In Doppler mode 
25 which is implemented by a Doppler processor 30, scanning is performed along a single 
scan line, and a spectrum of the blood velocity from a single spatial location is estimated 
from a substantially larger ensemble of data (typically, data from 64-5 12 pulses). Before 
displaying the processed image frame on a raster monitor or display 36, scan conversion 
is performed by a scan converter circuit 32, which converts the acquired ultrasound data 
30 from polar coordinates to the Cartesian coordinates used by the raster display. Post 
processing may optionally be applied by a post-processing circuit 34, to improve the 
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quality of the displayed image, as well as to combine the anatomy and flow images on 
the display. 

A tissue vibration processor 28 that is used to process the ultrasound data in the 
present invention is shown in FIGURE L The ultrasound data acquisition for tissue 
5 vibration processing and imaging is the same as that used for color-flow, so that 
quadrature data ensemble or color-flow data are input to the tissue vibration processor. 
However, instead of estimating blood flow velocity from these data, the tissue vibration 
processor estimates the instantaneous displacement of the scattering tissue from the phase 
of the complex received signal. This tissue motion is referred to as clutter in 
10 conventional color-flow imaging and is suppressed using clutter filters. Typically, 
cardiac pulsation, respiration and transducer motion each can contribute to an observed 
displacement or motion of tissue. Such motion is at a low frequency of a few Hertz or 
less. When internal bleeding is present, the tissue vibrates locally with a frequency 
ranging from a few tens of Hertz to a few hundred Hertz. By analyzing the frequencies 
15 of the different components of tissue motion, vibrations caused by bleeding may be 
distinguished from clutter caused by other sources of movement. The tissue vibration 
processor performs this analysis by decomposing the tissue motion into the dominant 
motion components and identifying any motion components that appear to be at a 
frequency higher than that of cardiac pulsation. 
20 It is contemplated that tissue vibration processor 28 can be implemented as an 

additional fixed-function circuit board or application specific integrated circuit (ASIC) in 
conventional ultrasound machines. Optionally, the tissue vibration processor can be 
combined with color-flow processor 26, since both process the same data ensemble. A 
standalone tissue vibration imaging device can be implemented with front-end 
25 subsystem 14, B-mode processor 24, tissue vibration processor 28, and scan 
converter 32. Those of ordinary skill in the art will appreciate that the tissue vibration 
processor can be implemented in software/hardware using one or more digital signal 
processors (DSPs) or alternatively, in an ASIC, or even on a conventional general 
purpose processor chip that access machine language instructions stored in a memory 
30 accessed by the processor to carry out the processing steps of the tissue vibration 
processor. 
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The computational power of ultrasound machines has increased significantly in 
recent years, benefiting from advances in processor technology. Thus, the additional 
computational burden of the tissue vibration imaging algorithms discussed below can be 
reasonably supported in modem ultrasound machines. Previously, a programmable 
5 ultrasound signal and image processing system suitable for use as the tissue vibration 
processor was developed using a new generation of high-performance multimedia 
processors to support all of the conventional processing modes, such as B, M, color-flow, 
and Doppler in software (Sikdar S, Shamdasani V, Gong L, Managuli R, Hayashi T, 
Mitake T, Kim Y. "A single mediaprocessor-based programmable ultrasound system," 
10 IEEE Trans Inf. Tech. Biomed 2003; 7:64-70), and subsequently, was shown to be useful 
in implementing tissue vibration processing in accord with the present invention. The 
main strength of a programmable system is the ease of development of new modes and 
applications such as tissue vibration imaging without the need for hardware 
modifications to conventional ultrasound machines. Integrated tissue vibration imaging 
15 using the software-programmable ultrasound system has thus been effectively and 
beneficially used for real-time visualization of vibrations in 2D ultrasound scans. 
Algorithms for Tissue Vibration Imaging 

In conventional color-flow imaging, the velocity of blood flow is estimated by 
computing the average phase difference between multiple ultrasound echoes (typically 
20 6-16 pulses) that are received from a sample volume. Echoes from moving tissue tend to 
have a significantly higher signal strength (typically 40 dB-60 dB higher), compared to 
the weak scattering from blood, and also have lower velocities. This high amplitude and 
low frequency tissue signal is commonly referred to as clutter and tends to bias the 
estimated blood flow velocity. Thus, clutter is suppressed using appropriate filters in 
25 conventional color flow imaging. The main components of clutter are cardiac pulsation, 
respiration, and transducer movement When blood flow eddies are present, any local 
tissue vibrations, e.g., those caused by the blood flow eddies at a bleeding site, will also 
be part of this clutter and would normally be suppressed in conventional ultrasound 
processing systems. 

30 The present invention separates the tissue vibrations from the remaining clutter 

and flow signals. In achieving this function, it was recognized that the tissue vibrations 
and clutter produce statistically independent signals that have different frequency content 
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While clutter due to cardiac pulsation and breathing typically occurs at 1 Hz or less, 
tissue vibrations typically occur at 50 Hz or more. Other noise sources are at 
substantially higher frequencies. Scattering from tissue is typically more coherent 
compared to the scattering from blood, because the tissue scatterers are more tightly 

5 bound together and tend to move as a group. Thus, compared to the clutter from other 
sources and tissue vibration signals, the blood flow signal typically has a much greater 
frequency bandwidth. Due to its weak signal strength and greater bandwidth, blood flow 
signals may be considered as noise compared to the stronger and more coherent tissue 
vibration signals for purposes of this invention. Therefore, tissue vibrations can be 

1 0 distinguished from clutter and blood flow based on spectral analysis. Spectral analysis of 
the phase of the received ultrasound echo can be used to separate the components of the 
scatterer motion, ignoring the scattered signal strength, whereas spectral analysis of the 
complex ultrasound echo considers both the signal strength and the motion components. 
Due to the limited number of temporal ultrasound samples (6-16 pulses) 

15 preferably used in the present invention, conventional Fourier transform and filtering 
approaches lack sufficient resolution for separating the signals. Therefore, this invention 
employs high-resolution spectral estimation techniques to carry out this function. Two 
high-resolution spectral estimation techniques were identified as suitable for this purpose, 
including eigendecomposition-based spectral estimation, which models the signal as an 

20 optimum set of orthogonal components, and autoregressive spectral estimation, which 
models the signal as the output of an autoregressive linear prediction filter driven by 
white Gaussian noise. Accordingly, three signal processing algorithms were developed 
for isolating tissue vibrations. The first algorithm is based on an 
eigendecomposition-based spectral analysis of the phase, of the received ultrasound echo; 

25 the second algorithm is based on an eigendecomposition-based spectral analysis of the 
complex ultrasound echo; and, the third algorithm is based on an autoregressive spectral 
analysis of the complex ultrasound echo. Since eigendecomposition is a 
computationally-intensive operation, an approximate eigendecomposition utilizing 
iterative QR factorization is used to develop a computationally-efficient algorithm. 

30 Signal Model 

To model the received signal from vibrating tissue, the tissue being imaged is 
approximated with S point scatterers having uniform motion and randomly distributed at 
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locations (F, = [r,]e r +[^,^ + s = l in a sample volume where 

{e r ,e v ,e^) denote the unit direction vectors in spherical coordinates. The instantaneous 
position of the scatterers,v(F,0 , is given by: 

F,(0 = \r, - d r (t)]e r + ty, - d v % + \4>, - 

5 where {d r {t),d v (0,^(0) denote the displacement as a function of time. If the scattering 
from the sample volume is uniform with a as the average scattering coefficient, then the 
scattering function of the sample volume is av(r,t) . The complex received signal from 
the sample volume, v(r,r), can then be modeled as a convolution of the pulse echo 
spatial impulse response, h pe <? ,t) , of a single point scatterer, the temporal response of 

10 the transducer, and the scattering function, av(F, t) . 

x(r) = x 0 (r)e J2 ' f ' T 

where the time indices r and t refer to "fast" time and "slow" time, respectively, f 0 is Ihe 
center frequency of the transducer, and «(r,r)is white thermal noise. Combining 
Eqs. (1), (2), and (3) results in: 

* C 

where c is the speed of sound and 2 fc"^(0) is the two-way pulse propagation time 

i c 

between the transducer and each point scatterer. The scatterer displacement for 
vibrations is small compared to the spatial size of the pulse echo spatial impulse response 
and the envelope of tire transducer response. Thus, the "slow" time variations in the first 
20 two terms may be neglected and Eq. (3) can be simplified to: 



(4) 



= A(r)e c +n(r,t) 



where A(r) is the complex amplitude of the scattered signal. It is apparent that the 
complex received signal is phase modulated with the instantaneous radial displacement. 
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If the scatterers in the sample volume are all vibrating radially in a simple 
harmonic fashion with peak displacement a 0 and frequency/**, the tissue displacement 
due to cardiac pulsation, breathing, and other tissue movement (i.e., clutter or noise) 
relative to the transducer is rf to (0- This motion will hereinafter be referred to as the 
5 "clutter motion." Then, the combined displacement can be considered to be a 
superposition, as follows: 

d r it) = d, iss (f) + a 0 sin(2^r vj ,r) (5 ) 

An ensemble of ultrasound pulses is transmitted in the same direction at a rate 
known as the pulse repetition frequency (PRF). Then, the complex received signal from 
10 the m th pulse transmission, y{r 9 m), is: 

y(T,m) = A(r)e J ° c ~ c +n(r,m) (6) 

where T PRF is the pulse repetition interval. The Fourier transform of the phase-modulated 
complex received signal is a Bessel series: 

Y(T,f) = A(j)\5(e J " c )*3{e ° ° ) 



= ^(r)|c(/)*^yV n (^(/-2^7' / , RF / vi6 )| + iV(r s /) (?) 
= A(r) "ZfJMctJ -27mT PRF f vib ) + N{r,f) 

n=~<x> 

15 where J, axe Bessel functions of the first kind, P = J0 0 , Sis the Dirac delta function, 

c 

and c(f) is the spectrum of the clutter motion (the clutter spectrum), and N(r,f) is the 
noise spectrum. 

FIGURE 2 illustrates a typical power spectrum 40 of the ultrasound signal when 
a tissue vibration is present. The spectrum includes multiple copies of the clutter 
20 spectrum separated by the vibration frequency, as indicated by Eq. (8), which is 
presented below. A low frequency peak 42 at/ to corresponds to the clutter spectrum, 
while symmetric peaks/, and/+ indicated respectively by reference numbers 44 and 46 
correspond to vibration, and puss, p+, and j?., are the corresponding peak powers. The 
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frequency peaks at f+ and f. are referred to herein as a "matching pair. 5 ' For 
small-amplitude vibrations, higher-order terms can be ignored; thus, most of the spectral 
energy will be present in the three frequency peaks, Z,^, and f+, respectively. Since, 

w P > the ratio of the power in the frequency peaks can provide an estimate of the 

vibration amplitude. Therefore, the vibration frequency and amplitude may be estimated 
from the power spectrum as follows: 



'power 
vib 



/,-/-[ 



<r=-r- F J J ^ t =-i (8) 



These estimators are referred to herein as the "spectral frequency estimator" and 
the "power ratio amplitude estimator." Alternatively, the vibration frequency and 
10 amplitude may be estimated from the residual phase {^(fc^^of the ultrasound signal 

after suppressing the effects of clutter motion. A coarse computationally-efficient 
estimate of the frequency of the dominant components, f vib , can be obtained by counting 
the zero crossings, Nzero, in the residual phase. This estimate can be further refined by 
interpolating the residual phase to compute the mean period of oscillation. The vibration 
15 amplitude may be estimated by the variance of the residual phase. These estimators are 
defined as follows: 



/* phase 
v 



9 vib 



a^r=^T var(^)) (9) 
W 0 



and are respectively referred to herein as the "zero-crossing frequency estimator" and the 
"phase variance amplitude estimator" 

20 For real-time tissue vibration imaging, only a short ensemble of ultrasound data 

(typically, 6-16 pulses or echoes) from each sample volume in a region of interest is 
available for processing. Conventional color-flow imaging systems utilize clutter 
filtering to suppress the clutter, while retaining the blood flow. However, due to the 
small number of temporal samples, the conventional clutter filtering-based approach, or a 

25 Fourier-based approach lacks sufficient resolution to discriminate between the tissue 
vibrations, blood flow, and clutter. A parametric approach that utilizes the characteristics 
of the vibration signal appears better suited to make this distinction. Three parametric 
approaches may be taken, based on the model of the ultrasound signal in Eqs. (6) and (7), 



WO 2005/037060 



PCTAJS2004/032427 



-16- 



including: (a) estimation of a pair of complex exponentials in noise; (b) autoregressive 
modeling; and, (c) decomposition of the phase of the ullrasound signal. In the following 
section, vibration detection algorithm based on these three parametric approaches are 
described in greater detail. One method of producing the complex ultrasound signal in 

5 Eq. (4) is a quadrature demodulation of the received ultrasound signal. An alternative 
method is to. compute the time delays producing the phase variations in Eq. (4) by 
processing the received RF ultrasound data using a cross correlation technique. 
Vibration Imaging Using Estimation of Co m plex Expo nentials in Noise 

Using the inverse Fourier transform of the Bessel expansion in Eq. (7), Eq. (6) is 

10 expanded, as follows: 



Thus, the ultrasound signal can be modeled as a sum of complex exponentials 
embedded in noise. As can be seen ftom the expression enclosed by square brackets in 
Eq. (10), vibrations correspond to matching pairs of complex exponentials. In contrast, 

15 the complex exponentials corresponding to clutter motion will typically not have such 
matching pairs of frequencies. The frequencies may be estimated using the 

root-MUSIC and ESPRIT algorithms (disclosed by P. Stoica and R. Moses in 
"Introduction to Spectral Analysis," Upper Saddle River, NJ: Prentice-Hall, 1997). 
Vibrations may then be detected using a matching peak criterion 

20 |/ + + /_ - 2/ <to | < F lhresll0ld , and the vibration amplitude and frequency can be estimated 
using Eq (9), which is set forth above. The steps of the algorithm are described in more 
detail below, in regard to FIGURE 3A. Based on this criterion, vibrations can be 
detected and distinguished from clutter motion. Any blood flow signals may be 
considered as part of the noise spectrum 

25 FIGURE 3 A illustrates a flow chart 50 that shows the logical steps involved in a 

first algorithm for estimating the tissue vibrations based upon a pair of complex 
exponentials in clutter or noise, that are normally excluded from color-flow processing. 
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The procedure begins with a quadrature-demodulated ensemble of 2D ultrasound 
data 52. 

In a step 54, a correlation matrix is estimated from the color-flow data. In a 
step 56, the correlation matrix is employed to carry out an eigendecomposition, 
5 producing a signal subspace 58 and a noise subspace 60. Using the signal subspace and 
the noise subspace, the frequency of the dominant components is estimated in a step 62, 
by employing the root-MUSIC and ESPRIT algorithms, as noted above. A decision 
step 64 then determines if each dominant component is double-sided, while a step 66 
estimates the vibration amplitude and frequency of each dominant component using 

10 Eq. (8). If a dominant component is not double-sided (i.e., is not a potential tissue 
vibration component), decision step 64 returns a "zero," while if the dominant 
component is double-sided, the decision step return a one. A multiplier 68 then 
multiplies the output of decision step 64 times the vibration amplitude and frequency 
estimate for the dominant component, yielding a null if the dominant component is not a 

15 tissue vibration component and returning the estimate of vibration amplitude and 
frequency of the dominant component otherwise. A median filter 70 then filters isolated 
falsely-detected vibrations and other undesired noise from the results, so that the 
remaining vibration image indicating a bleeding site is displayed in a step 72. 
Vibration Ima ging Using an Autoregressive Signal Model 

20 The ultrasound signal from vibrations can be modeled as the output of ap^-order 

autoregressive linear prediction filter with white Gaussian noise having a variance o^, as 
the input, as follows: 

y(r, ™) » Yi a >n~* (*)y(?> »»-*) + n(r, m) 

(11) 

The linear prediction coefficients, a*(i), can be computed using either a 
25 least-squares minimiz ation of the prediction errors or using the computationally-efficient 
Burg algorithm, as explained by Stoica and Moses in the above-referenced paper. A 
high-resolution spectral estimate can then be obtained from this autoregressive model as 
follows: 
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From Eq. (8), the presence of symmetric matching pairs of frequency peaks in the 
power spectrum around the clutter motion peak may be detected as a vibration. As 
before, any flow signals may be regarded as noise. For ultrasound signals from 
vibrations, the power spectrum in Eq. (13) will have frequency peaks 0^,/+,/) at the 



1 + 



Vibrations can be 



local minima of the polynomial A(r, f) = 

detected using a matching peak criterion |/ + + /_ - 2/^ | < F 0tres}wld , and the vibration 
amplitude and frequency can be estimated using Eq. (9). The steps of this algorithm are 
described in more detail below, in connection with FIGURE 3B. 

10 As shown in a flow chart 80 in FIGURE 3B, the second alternative algorithm 

also begins with quadrature-demodulated ensemble data set 52. In a step 82, reflection 
coefficients are computed for each ensemble of the quadrature-demodulated data. Using 
the reflection coefficients, linear prediction coefficients are determined in a step 84. In a 
step 86, the power spectrum is estimated from the linear prediction coefficients and the 

15 peaks in the power spectrum are detected. Again, decision step 64 determines if the 
peaks thus identified are for tissue vibration by determining if they are double-sided and 
returning a zero if not, and a one, if so. Also, step 66 provides for estimating the 
vibration amplitude and frequency at each of these peaks, and the results from decision 
step 64 are multiplied the estimated amplitude and frequency in multiplier 68. Median 

20 filter 70 is then applied to the results and the filtered image data are displayed as a 
vibration image, in step 72. 
Vibration Imaging Based on Phase Decomposition 

Another model of vibrations can be based on the phase ultrasound signal. As 
shown in Eq. (7), vibrations will produce an oscillatory signature in the phase, which will 

25 typically not be present in the case of clutter motion. Although flow signals may have 
oscillatory phase, the echoes from vibrating tissue are expected to be more coherent than 
those from flow. Thus, their phase may be modeled by a smaller number of dominant 
components. Accordingly, a vibration detection algorithm can also be based on 
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decomposition of the phase of the ultrasound signal into its do min a nt components and 
testing for oscillatory phase. Alternatively, instead of using quadrature-demodulated 
ultrasound data, the phase can he estimated from RF ultrasound data by estimating the 
time delays between a pair of RF ultrasound data, 

5 Any linear time-varying motion is first suppressed by down mixing the ensemble 

of 2D ultrasound data with the mean clutter velocity, estimated using the conventional 
autocorrelation method. The phase of the ensemble of 2D ultrasound data is then 
computed, and the mean phase is subtracted to suppress the effect of the stationary echo. 
The residual phase is then decomposed into its dominant components using a method 

10 similar to principal component analysis. The first step of the decomposition involves the 
estimation of the correlation matrix of the residual phase using the modified covariance 
method (Maiple, 1987). An approximate eigendecomposition can then be performed 
vising iterative QR factorization of the correlation matrix. The approximate eigenvalues, 
JLi 9 may be estimated by the diagonal elements of the upper triangular matrix Rk after the 

15 A* iteration. The eigenvectors are arranged in order of decreasing eigenvalues. The 
eigenvalues are a measure of the signal energy contributed by the corresponding 
eigenvector. Thus, the fraction of the total signal energy contained in the p dominant 

components can be estimated using E p = ^ — . Therefore, noise and blood flow can 

i=l 

be suppressed using a threshold criterion, E p >E t hre S hoid* To further separate tissue 
20 vibrations from clutter motion, the fact that tissue vibrations have a higher frequency 
compared to clutter motion is applied Vibrations can then be separated from clutter 
using a frequency threshold criterion f vib > F lhreshold , where F threshold is chosen so that at 
least one half of one period of the vibration is contained in an ensemble. The vibration 
frequency and amplitude may be estimated using Eq. (10). The steps of this algorithm 

25 are described in more detail below, in connection with FIGURE 3C 

A flow chart 90 illustrates the logical steps of the third algorithm. Again, starting 
with the quadrature ensemble data, a step 92 down mixes the quadrature data with a 
mean clutter velocity, which is determined from the data. A step 94 provides for 
computing an unwrapped phase of the quadrature ensemble or color flow data, and then 

30 subtracting the mean clutter velocity from the unwrapped phase, resulting in a residual 
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phase. As explained above, instead of determining the phase from quadrature- 
demodulated data, the phase can be determined from RF ultrasound data by estimating 
time delays between a pair of RF ultrasound data. Using the residual phase, a step 96 
estimates a correlation matrix, which is then used to carry out a QR factorization in a 

5 step 98, yielding an eigenvalue estimate 100, an eigenvector estimate 102, and a vibration 
amplitude and frequency estimate 104, which are determined using Eq. (10), as noted 
above. Using the eigenvalue estimate, a decision step 106 determines if the total energy 
contained in the p dominant component is greater that a predefined threshold, T. If so, 
decision step 106 returns a zero if not, and a one if so. Similarly, a decision step 108 

10 determines if the estimate eigenvector has a frequency that is greater than a predefined 
threshold, F. If so, decision step 108 returns a one, and if not, a zero. The results of 
decision steps 106 and 108, and the estimated vibration amplitude and frequency of the 
dominant components are then multiplied together by a multiplier 1 10, so that if either of 
the decision blocks has returned a zero, the result is null, but if neither has returned a 

1 5 zero, the estimated vibration amplitude and frequency from step 104 are returned. Again, 
median filter 70 is applied to the estimated amplitude and frequency, providing filtered 
results that are displayed as the vibration image, indicating a bleeding site, in step 72. 
Detectable Vibration Amplitudes and Frequencies 

hi experiments using a physical phantom model, tissue vibrations with a peak 

20 amplitude of about 1 |wn have been accurately detected. The minimum detectable 
vibration amplitude depends upon the noise level and dynamic range of the phase of the 
received ultrasound echo. In modem ultrasound machines, the phase can have a dynamic 
range of 96 dB or more (for 16-bit quadrature-demodulated data) and the signal exceeds 
the electronic and thermal noise level by typically 80 dB or more. Therefore, from 

25 Eq. (4), vibrations as small as 50 nm may theoretically be detected using ultrasound. 
Practically, the attenuation of the ultrasound signal will reduce the dynamic range and 
limit the minimum detectable amplitude in deep tissue to -0.5 jam. 

The detectable vibration frequencies depend upon the choice of PRF, i.e, on Fprf* 
A PRF that is too low compared to the vibration frequency would lead to aliasing, while 

30 selecting a PRF that is too high will fail to detect low-frequency vibrations. A vibration 
can be detected only if at least half of one vibration cycle is captured within the temporal 
window corresponding to an ensemble. Thus, all vibrations with frequency between 
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F F 

— ^ and p ^ can be detected theoretically without aliasing for an ensemble size E. 

Since vibrations can be broadband, a high-frequency vibration interrogated at a low PRF 

value can be mistaken for noise using this algorithm. Thus, for better sensitivity, it is 

desirable to select a PRF and an ensemble size so that only a few periods of the vibration 

kF 

5 are included in an ensemble. Accordingly, the maximum detectable frequency is — 

E 

when k periods of the vibration are included in an ensemble. A simulation and phantom 
experiments that were carried out indicate that reliable detection may be performed using 
only one half to six vibration periods during the interrogation period. For example, with 
a PRF of 1 kHz and an ensemble size of 1 6, vibrations with frequency between 3 1 .3 Hz 

1 0 and 375 Hz may be reliably detected. 
Determination of Bleeding Rate 

Since the tissue vibrations are produced by the blood flow eddies, the frequency 
of the tissue vibrations is the same as the frequency of the eddies. The Strouhal number 
(S) relates the frequency of eddies (f vi b) produced by an orifice to the diameter of the 

1 5 orifice (D\ and the flow rate through the orifice (£7), according to: 

S= fvib * D (13) 
U 

Since the Strouhal number remains constant, the vibration frequency can be related to the 
diameter of the orifice and the flow rate through the orifice. In addition, the energy in the 
eddies (E) and thus, the amplitude of the tissue vibrations (a v &) 5 is directly proportional to 
20 the flow rate, as follows: 

Eoza 2 vib ocU 2 (14) 
KxD 2 xU 

The bleeding rate Q = can thus be determined from the frequency of the 

vibrations and the amplitude of the tissue vibrations. 
Sources of Artifacts 

25 In color-flow data acquisition, interrogation along each scan line is performed for 

only a brief period of time. Vibrations are transient, with typical durations of 
10 ms-100 ms. Thus, there is a possibility that some vibrations may not be interrogated. 
Since the vibrations typically have a relatively large spatial extent and repeat every 
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cardiac cycle, it is unlikely that the vibrations will be missed entirely; however, the 
spatial extent of the vibrations visible in the image may be only a part of the true spatial 
extent. By appropriately choosing the PRF and the region of interest, such discrepancies 
may be minimized. 

5 Other artifacts may be falsely detected as vibrations. Transducer motion may 

introduce additional frequency peaks in the clutter spectrum and may cause false 
detections; however, these false detections may be minimized by using a trained 
sonographer to perform the sc anning . Vibrations in the tensed skeletal muscle of the 
sonographer, and any ambient vibrations may be detected in the vibration image. In 

10 addition, the high-resolution spectral estimation methods may produce spurious peaks 
that can be falsely detected as vibrations. Such artifacts can be easily distinguished from 
pathological vibrations, which are expected to be correlated with the anatomy and 
periodic with every cardiac cycle. These artifacts can be also avoided if additional 
temporal samples are available. Any vibrations displayed in the vibration image should 

1 5 therefore be confirmed with the vibration spectrum by placing a Doppler sample volume 
at the location of the peak intensity. 
Comparison of the Algorithms Derived from Modeling 

The ability of the proposed algorithms to detect vibrations was evaluated using a 
simulation model. Simulations show that subspace-based algorithms such as MUSIC 

20 and ESPRIT have high sensitivity (96%) and specificity (98%) for detecting narrowband 
vibrations in the presence of clutter as well as blood flow and are robust even when 
broadband vibrations are present. For narrowband vibrations, an algorithm based on an 
autoregressive model has a slightly improved specificity (99%), a comparable sensitivity, 
and is robust to broadband vibrations. The phase decomposition-based algorithm has a 

25 slightly lower sensitivity (93%) and specificity (98%), but is more robust to broadband 
vibrations. 

The computational requirements of the proposed algorithms are shown below in 
Table 2. The subspace-based algorithms (MUSIC/ESPRIT) have a computational 
requirement that is highly dependent on the choice of model order. In these algorithms, 
30 eigendecomposition is the most computationally-intensive task. The 
autoregression-based algorithm is less computationally intensive and the computational 
requirement is less dependent on the model order. In this case, the computation of the 
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FFT for spectral estimation is the most computationally-intensive task. The phase 
decomposition method is the least computationally intensive, since it involves operations 
on real signals only. Thus, the phase-decomposition algorithm is most suitable for 
real-time implementation. 

Table 1 - Field II Simulation Parameters 



Center frequency 
PRF 

Transducer excitation 
Transducer impulse response 
/number 

Number of elements 
Transducer height 
Element pitch 
Element kerf 
Transmit aperture 
Receive aperture 
Transmit focus 
Receive focus 
Elevation focus 
Mathematical element size 
Sampling frequency 
Sound velocity 
Number of scan lines 
Number of ensembles 



5MHz 

500 Horizontal 
5-period sinusoid 

Hanning-weighted 2-period sinusoid 
2 

192 

15 mm 
0.4 mm 
0.03 mm 
25.6 mm 
25.6 mm 
40 mm 

30 mm to 100 mm in steps of 10 mm 
20 mm 

0.37mm x 1.5 mm 

105 MHz 

1540 m/s 

32 

10 



Table 2 - Computational requirement (million operations/sec) for real-time 



imaging at 10 frames/s with 32 scan lines, 256 samples/scan line, and ensemble 10 
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Algorithm 


Model Order ~] 






P=4 


MUSIC 




3631 


7653 


ESPRIT 




2218 


6455 


AR 




1606 


1630 


Phase-decomposition 


181 


489 


1107 


Color flow 


89 
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In Vivo Empirical Tests of Arterial and Organ Bleeding in A nimal Model 

FIGURE 5A illustrates a color-flow image of a punctured femoral artery. A 
dotted line box indicates the region of interest An arrow indicates the approximate 
location of a puncture. In FIGURE 5B, a vibration amplitude image of the punctured 

5 femoral artery is illustrated Although not visible in color in this grayscale image, the 
vibrations at the puncture site in the original color image appear green relative to the 
shades of gray in the surrounding tissue. Similarly, in FIGURE 5C, a vibration 
frequency image shows low-frequency vibrations in color at the puncture site, relative to 
shades of gray for the surrounding tissue. 

10 FIGURE 4 is a schematic view 120 that illustrates how ultrasound transducer 12 

is employed to produce color-flow data for an incised organ 124. To couple acoustic 
signals from and to the ultrasound transducer, a gel standoff 122 is disposed between the 
ultrasound transducer and the exposed organ. The gel standoff provides an efficient 
couple between the tissue of the organ and the transducer. In this manner, tissue 

1 5 vibrations caused by bleeding through an incision 126 can readily be produced in accord 
with the present invention. 

In FIGURE 6A, a color-flow image of an incised spleen is shown. An arrow 
indicates the pooled blood from an incision, where the blood pool has accumulated 
between the organ and the gel stand-off. Some blood flow is observed at the bleeding 

20 site as indicated by an arrow. FIGURE 6B shows the vibration amplitude image, where 
vibrations surrounding the bleeding site are indicated by an arrow. A vibration frequency 
image in FIGURE 6C shows the low-frequency vibrations in color to indicate the 
bleeding site. 

FIGURE 7A shows the color-flow image of an incised liver in a rabbit. The liver 
25 was surgically exposed and a triangular incision was made with a scalpel, causing 
moderate bleeding. The liver was imaged though a gel standoff. The location of the 
incision is clearly seen in the B-mode image (not shown), and some flow is observed in 
the artery feeding the bleed. However, from the color-flow image, it is not possible to 
identify the bleeding site at the incision. The vibration amplitude image is shown in 
30 FIGURE 7B. Vibrations are clearly observed surrounding the incision from which 
bleeding is occurring, as indicated by an arrow. The vibration frequency image in 
FIGURE 7C shows low-frequency vibrations, and the bleeding site is again clearly 
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visible (although not so apparent in this grayscale image as in the original color images). 
The vibration image and color-flow images therefore provide complementary 
information and a clear correlation between the spatial location of the vibration, and the 
underlying anatomy can be observed Thus, a combination of B-mode, color-flow, and 
5 tissue vibration images can be a powerful diagnostic tool for detecting the location of 
internal bleeding. 
Summary 

Empirical evidence demonstrates the feasibility of real-time ultrasound imaging 
of low-intensity local vibrations in the vessel wall and surrounding tissue associated with 

10 punctured arteries and organ bleeds. Several algorithms based on parametric signal 
decomposition and spectral estimation have been developed for ima g ing small-amplitude 
tissue vibrations using as few as 10 temporal samples. The vibration amplitude and 
frequency can be estimated accurately, and real-time tissue vibration imaging has been 
implemented on an ultrasound machine with a software-programmable subsystem. 

15 Tissue vibration imaging can provide additional diagnostic information that is 

currently not available to the clinician using conventional tools. An ultrasound device 
with tissue vibration imaging capability can become a useful screening and diagnostic 
tool for paramedics and trauma centers for rapid diagnosis and localization of active 
internal bleeding. Such real-time localization of bleeding can then be employed for 

20 targeting and evaluating the effects of HIFU and other therapies intended to stop the 
bleeding. 

Although the present invention has been described in connection with the 
preferred form of practicing it and modifications thereto, those of ordinary skill in the art 
will understand that many other modifications can be made to the present invention 
25 within the scope of the claims that follow. Accordingly, it is not intended that the scope 
of the invention in any way be limited by the above description, but instead be 
determined entirely by reference to the claims that follow. 



